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We analyze the stability of two charged conducting spheres orbiting each other. Due to charge 
polarization, the electrostatic force between the two spheres deviates significantly from as they 
come close to each other. As a consequence, there exists a critical angular momentum, Lc, with a 
corresponding critical radius rc- For L > Lc two circular orbits are possible: one at r > rc that is 
stable and the other at r < rc that is unstable. This critical behavior is analyzed as a function of 
the charge and the size ratios of the two spheres. 



I. INTRODUCTION 

Although Coulomb's law mimics Newton's law of grav- 
itation in its behavior, no electrostatic orbits are 
found in nature while gravitational orbits are abundant. 
The reason is that electrostatic forces are much stronger 
than gravitational forces and form bound systems at 
much smaller, atomic length scales. Electrons orbit the 
nucleus in Bohr's atomic model, but quantum mechanics 
shows that atomic orbitals are much more complex. 

Creating macroscopic electrostatic orbits is an exper- 
imental challenge because of the requirements of high 
voltages and nearly zero gravity. In addition, there are 
practical difficulties such as leakage of charge from the 
spheres to ambient air and sparking to nearby grounded 
objects. In August 2006 Rhodes College undergraduates 
successfully demonstrated electrostatic orbits^ by mak- 
ing a small graphite coated negatively charged Styrofoam 
sphere orbit a larger stationary aluminum sphere in con- 
ditions of microgravity created aboard the NASA Mc- 
Donnell Douglas C9B aircraft. In July 2008 Rhodes Col- 
lege undergraduates created binary orbits between two 
oppositely charged equal sized graphite coated Styrofoam 
spheres.'^ 

In these experiments the spheres crash into each other 
if they come too close. The reason is the redistribution 
of charge on the two spheres due to polarization which 
causes the force between the two spheres to deviate signif- 
icantly at shorter distances from 1 /r^ behavior. At some 
separation r the resultant force is able to overcome the 
"centrifugal repulsion" due to the angular momentum of 
the orbiting objects thus making the orbits unstable. In 
this paper we analyze the stability of electrostatic orbits 
taking polarization effects into account. 

The electrostatic force between two conducting spheres 
with arbitrary charges and radii is a fundamental prob- 
lem that has generated considerable interesti^ Classic 
books on electromagnetism by Maxwell,'* Jeans?^ and 
Smythe^ include a detailed analysis of this problem. 
Maxwell credits Poisson'*-- as the first to solve the prob- 
lem and Thompson (Lord Kelvin) for a "greatly simpli- 
fied" treatment lii^ More recently, Soules^ did a numeri- 
cal analysis of the force between two oppositely charged 



spheres of unequal sizes to show that polarization effects 
dominate at short distances. 

The main idea behind the solution to this problem is 
the method of images. A complete analysis of the prob- 
lem requires an infinite number of image charges. How- 
ever, as shown numerically in Refs.fs! and'o', only a finite 
number of image charges are needed to compute the force 
to a desired degree of accuracy. This number increases 
rapidly as the spheres come closer. For the analysis in 
this paper we take 120 image charges on each sphere into 
account. 

Because the interaction between the two spheres de- 
pends only on their center-to-center distance r, we can 
incorporate the angular motion of the objects into a 
centrifugal term, and reduce the problem to a one- 
dimensional problemiiSiii For interaction forces of the 
type F{r) ~ l/?"", stable orbits are possible for n < 3. 
For n > 3 the repulsion is unable to overcome the 
1/r" attraction at short distances; at large distances the 
repulsion becomes stronger than the attraction, and thus 
the stability is lost at both ends. 

For our problem stability is not an issue at large dis- 
tances because the Coulomb force dominates the polar- 
ization effects for large r. At small r polarization ef- 
fects become significant and higher and higher powers of 
1 /r are needed to describe the interaction accurately. In 
Sec. [n] we analyze the stability of a point charge orbiting 
a charged conducting sphere, and in Sec. lIIII we study the 
more general case of two spheres with arbitrary charges 
and radii. We ignore energy loss due to radiation from 
the centripetal acceleration of the orbiting charge. By us- 
ing Larmor's formula^^ it can be shown that the lifetime 
of the system due to radiation loss will be much larger 
than its orbital periodj^^ 



II. POINT CHARGE ORBITING A 
CONDUCTING SPHERE OF CHARGE 

Consider a conducting sphere A with charge Q and 
radius a with a point charge —q of mass m orbiting 
around it. Without loss of generality we let A be fixed 
and centered at the origin. Binary orbits between two 
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bodies about their center of mass can be reduced to an 
equivalent one body problem by switching to relative 
coordinates.— This problem requires only one image 
charge; it is exactly soluble and helps to understand the 
qualitative features of the more general case discussed in 

sccnni 

Let the distance of —q from the origin be r. Due to 
the presence of —q, the charge on A is redistributed and 
can be considered as being composed of a charge Qo a-t 
xq = and Qi — aq/r at xi = a^/r, where xq and xi are 
measured from the origin along the line joining the two 
charges. Because Q is fixed, we have Qq — Q ~ aq/r. 

If we use the principle of superposition, the net force 
on the point charge can be written as 
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where k — I/Attco. Integrating Eq. ([T]) with respect to r 
(from oo to r) gives the electrostatic interaction energy 
of the system as 
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The effective potential energy for the system is 
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where L is the total angular momentum of the system, 
and the first term in the energy is the centrifugal term. 
For numerical and analytical calculations it is convenient 
to define the dimensionless variables 
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We first set Q = \ and plot C/cff(r) for several values 
of L^. For P = 2.83 (see curve (a) in Fig. [J), the ef- 
fective potential has two extremum points correspond- 
ing to two circular orbits: the one at f = 2.32 is sta- 
ble {(PtJes/df"^ > 0), and the one at f = 1.87 is un- 
stable {d^UcS / df^ < 0). This smaller unstable orbit is 
purely due to polarization effects because at sufficiently 
close distances the polarization effects are so strong that 
they overcome the centrifugal repulsion. As we decrease 
th^ value of we find that there exists a critical value 
Lc = 2.7715 where the two extrema merge into a single 
point at the critical radius, fc = 2.0704 (see curve (b) in 
Fig. [IJ. There are no stable or unstable orbits below this 
critical value of (see curve (c) in Fig.[T]). 
To obtain the critical parameters we set 
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and 
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FIG. 1: !7cff(f) versus f for (a) = 2.83. The arrows in- 
dicate stable and unstable circular orbits at f = 2.32 and 
f — 1.87 respectively, (b) The critical angular momentum 
Zc = 2.7715. The corresponding critical radius (indicated by 
the vertical dashed line) is fc = 2.0704. (c) = 2.70 with 
no possible orbits. 
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By eliminating from these conditions, we can write 
the equation for as 
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The substitution of IJe (defined through Eqs. ([2]) and 
(gl)) into Eq. jH) gives 
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On the left-hand side of Eq. ([7]) is the Coulomb force due 
to the first term in Ue from Eq. and on the right- 
hand side is the expression for the polarization force.— 
The point at which these two are equal determines the 
critical radius. Note that the polarization force is inde- 
pendent of Q because the charge on the sphere is unable 
to polarize the point charge. It goes to oo in the limit 
f — !■ 1 as expected. For large f it goes as 1/f^ because 
the orbiting point charge creates a dipole of order l/f^ 
in sphere A that interacts with the point charge with its 
dipolar field that goes as 1/f^. 

The solution of Eq. ([7]) for fc gives 



1 + 
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where z = {2Q) . In the limit Q ^ oo Eq. ([5]) gives 
fc = 1 (or Tc = a) because the Coulomb force in Eq. ([7]) 
dominates the polarization force for f > 1. In the other 
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limit Q — > 0, Tc — > oo as XjQ^I^ . The polarization force 
dominates the Coulomb force for f ^ in this limit, 
and Eq. ([7|) in the leading order becomes Qjr^ ~ 
Solving for fc gives the XjQ^I'^ behavior. 

The substitution of in Eq. (|5ap or Eq. (j5b[) gives 
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In the limit Q ^ co, L\ approaches oo as the first power 
of Q. This result can be understood by setting the cen- 
trifugal repulsion L\lf\ equal to the Coulomb force^^ and 
noting that = 1. For Q ^ 0, approaches as Q^^^ . 
This result follows from the same argument as before and 
noting that ~ XjQ^^^ in this limit. 



III. ORBITING CHARGED CONDUCTING 
SPHERES OF UNEQUAL SIZE 

We now consider the more general case of two conduct- 
ing spheres A and B with charges qi = Q and q2 = —q 
and radii a and b, respectively. To calculate the effec- 
tive potential energy in Eq. ([3|) we need the electro- 
static interaction energy between the two spheres. As 
noted by Maxwell, this calculation requires "an intricate 
investigation"— because the two spheres create an infinite 
series of image charges on themselves and on each other. 

The total electrostatic energy of A and B (including 
their self energies) is given hy^ 

Wv{r) = liciiV,^ + 2c,2ViV2 + C22Vi), (10) 

where Vi and V2 are the voltages of A and B respec- 
tively, and the coefficients Cij are their self and mutual 
capacitances defined by the relations 
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qi = ciiVi -I- C12V2 

92 = C21V1 -I- C22V2. 



We use these relations to express the energy in Eq. JTO] 
in terms of the charges of the spheres as 



^?(^) = + 2pi2qiq2 



P22ql) 



(12) 



where the matrix pij is the inverse of the 2x2 capacitance 
matrix Cy. 

To calculate the coefficients Cij the interaction between 
the spheres is ignored as a zeroth order approximation, 
giving a charge of aVi/k on A and bV2/k on B. We 
then use the method of images to calculate successive 
corrections to the charges on each sphere so that they 
remain fixed at their given potentials Vi and V2. The 
infinite series of charges on the spheres is summed and 
set equal to their net charges qi and (72 and the coeffi- 
cients of capacitance are identified from Eq. ([TT|) . Pois- 



son and Kirchofi^ calculated these coefficients and ex- 
pressed them as 
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^2 (a + ba)^ , 2 (b + aa)^ 

^ = ^2 ^'^d ^ = ^2 ■■ 

and a is the smaller root of the quadratic equation 
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In the limit r — > 00, cn — a/k, C22 — b/k^ and C12 = 
C21 = as expected. We truncate the series for Cy in 
Eq. (fT5|) at n = 60. This truncation takes into account 
the contribution of 120 image charges on each sphere in 
addition to the charges at their centers. 

From the capacitance matrix Cij we calculate its inverse 
matrix and substitute it into Eq. and expand 
in powers of l/r to obtain the total electrostatic energy 
(divided by A:) as 

Wq{r) ^ gf _ <M2 _ fc^gi + a^g| _ b^ql + a^gj 

k 2a 2b r 2r^ 2r^ 



r 2r^ 
2a%^qiq2 _ b'^qf + a' qj 
7^ 2r8 
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After subtracting the self energies q\l2a and g|/26 of A 
and i?, and then setting qi = Q and g2 = —q, we expand 
their electrostatic interaction energy (divided by kq'^/a) 
in powers of 1 /f to obtain 
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where s = b/a is the size ratio of the spheres. This 
form of the interaction energy lets us easily calculate the 
derivatives and double derivatives needed in the follow- 
ing. Including 120 image charges on each sphere allows 
us to carry out this expansion to order (l/f)^2i and gives 
numerical values for the force {—(IUe I df) accurate to 
e < 10"*^ for r > l.l(a + b) and 0.01 < Q < 100. The 
quantity e is the fractional change in the force due to the 
inclusion of one additional term in the approximation of 



Cij in Eq. p3| . 

The rest of the treatment mirrors that in Sec. nil With- 
out loss of generality we assume s < 1 and let A be fixed 
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at the origin and B be the orbiting sphere. The effective 
potential energy is defined using Eqs. ^ and (g]). To 
find the critical radius we substitute IJe from Eq. p7|) 
into Eq. ^ and obtain 



4(s3g2 + l) 12(s 
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Equation (fT8| is the general form of Eq. ([7]) and takes 
into account the finite size of the orbiting sphere. At 
the critical radius the Coulomb force on the left has to 
balance the infinite series for the polarization force on the 
right. We solve Eq. (fT8|) numerically for fc as a function 
of the charge ratio Q for several values of the size ratio 
s. The results are plotted in Fig.O 
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TABLE I: The values of the charge ratio, Qmin, at which the 
minimum critical radius occurs are listed (see second column) 
versus the size ratio s. The third column gives the value of 
this minimum critical radius and the fourth column gives the 
corresponding value of For smaller size ratios the fixed 

sphere can hold a greater charge than the orbiting sphere, and 
the hmiting orbits are smaller and faster as indicated by the 
last two columns. 
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FIG. 2: The critical radius, fc, plotted against the charge 
ratio Q for size ratios of (a) s = 1, (b) s — 3/4, (c) s = 1/2, 
(d) s = 1/4, and (e) s = 0. All curves with non-zero size 
ratio asymptotically converge to the s — curve which goes 
as in the limit Q — > 0. In the limit Q ^ oo they go 

as sQ^^^ except for the s = case which goes asymptotically 
to 1. 

For Q ^ Eq. ^TE\i becomes independent of s, and 
we obtain the s = case discussed in Sec. |lTl In this 
limit all the fc versus Q curves asymptotically approach 
the s = curve, which goes to oo as For Q 

oo, due to the terms in Eq. (fTSl) . the polarization 
force dominates the Coulomb force for any finite value 
of f. The equation for the critical radius in this limit is 
Q/f^ — ■^s^Q^ /f^- Thus, fc goes to oo as sQ^^^ as seen 
in Fig. O Note that for all s ^ 0, the critical radius has 
a minimum, fc,min, with respect to Q. This minimum is 
expected because fc goes to oo at both ends of the (0, oo) 
domain for Q. The numerical values of fc,min are listed 
in Table [J along with the corresponding values of 
and the charge ratio Qmin at which they occur. 

For s = 1 the minimum critical radius is at Qmin = 1 as 
expected, because the result has to be symmetric with re- 



spect to the inversion Q 1/Q {q ^ Q). Equation (fT5)) 
is invariant under this transformation as be seen by di- 
viding by Q on both sides. For s < 1 the Q ^ 1/Q sym- 
metry is broken. As s becomes smaller, the minimum 
critical radius decreases and the value of Qmin at which 
it occurs increases (see TableUand Fig. [2]). Qualitatively, 
as the size ratio decreases, the larger fixed sphere is able 
to hold a greater charge without excessively polarizing 
the smaller orbiting sphere. This leads to smaller stable 
orbits with higher angular momentum. 

The numerical solution of fc can be substituted into 
Eqs. ([5a|l or (|5bp to obtain the corresponding solution of 
L^. The results are plotted in Fig. [3] In the limit Q — > 0, 
Eq. (fTS)) is independent of s, and all the curves approach 
zero as Q^^^, which is the same behavior as found for 
a point charge. The Q — > oo limit can be obtained by 
setting the centrifugal force L'^/fl equal to the Coulomb 
force Q/fc- In this limit fc ^ sQ^/^ and thus ^ oo 
as sQ4/3. 

It it interesting to compare this theory to the successful 
orbits reported in Table I of Ref. [l]. The size ratio of the 
spheres for these experiments was s = 1.6 cm/6.5 cm = 
0.25, and the charge ratio, Q, (estimated from the given 
voltageai^) ranged from 5.0 to 5.8. For these values the 
critical parameter is between 10.3 and 11.8 and fc 
is 1.7. The successful orbits had initial values— of 13 < 
< 42 and 2.2 < f < 2.8, which are above the threshold 
needed for stability. 



IV. CONCLUSIONS 

We analyzed the stability of orbits between two 
charged conducting spheres due to their electrostatic at- 
traction. For a two body gravitational system only one 
stable circular orbit is possible for a given value of angu- 
lar momentum. In contrast, for an electrostatic system 
two circular orbits are possible for a fixed angular mo- 



5 




FIG. 3: The critical angular momentum, L^i plotted against 
the charge ratio Q for size ratios of (a) s = 1, (b) s — 3/4, 
(c) s = 1/2, (d) s = 1/4, and (e) s = 0. All the curves 
approach zero as Q^^^ independent of the size ratio. In the 
limit Q oo they go as sQ^^^ except for the s = case which 
is linear in Q in that limit. 



possible stable and largest possible unstable orbits. The 
critical parameters in dimensionless form depend only on 
the charge and the size ratio of the two spheres. 

Analyzing the critical parameters with respect to the 
charge and size ratios yields a minimum critical radius 
for every non-zero value of the size ratio. This minimum 
critical radius occurs at the charge ratio of unity when 
the spheres are equal in size. As the size ratio decreases, 
the charge ratio where the minimum critical radius occurs 
increases while the value of the minimum critical radius 
decreases. In the limit of zero size ratio (point charge), 
the charge ratio approaches c» and the minimum critical 
radius approaches the radius of the stationary sphere. 

We suggest two possible undergraduate research 
projects as follow ups to the work in this paper. Al- 
though zero gravity conditions are ideal for observing the 
stability of electrostatic orbits, ground based experiments 
can be performed by hanging the orbiting sphere from a 
high ceiling. The orbits can also be simulated using the 
equations of motion and the full electrostatic interaction 
(including polarization effects) between the two spheres. 



mentum. The larger of these two orbits is stable and the 
smaller orbit is unstable. 

We also found that there exists a critical angular mo- 
mentum below which no stable or unstable orbits are 
possible. The critical radius corresponding to this critical 
angular momentum is the limiting value of the smallest 
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The ratio of the lifetime to the orbital period is of order 

c^/v^, where v is the orbital speed and c is the speed of 

fight. 

The polarization force, defined as the right hand side of 
Eq. ([7}, comes from substituting the non-Coulombic terms 
in Eq. ((2]) into Eq. (|6]). It is not equal to the total force 
minus the Coulomb force. 

At the critical point the centrifugal force is of the same 
order as but is not exactly equal to the Coulomb force. 
Ground based tests showed that the orbiting sphere re- 
ceived only about 80% of the maximum possible charge at 
that voltage due to the small but finite size of the electrode 
which "shared" the charge with the sphere. 
^"^ In addition to the parameters listed in Table I of Ref. [l| 
the launch velocities obtained from the video of the orbits 
were used in calculating the values of . A link to the 
video is provided in Ref. [J- 



